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1 Introduction 

The coupling of a system to its environment is a recurrent subject in this collec- 
tion of lecture notes. The consequences of such a coupling are threefold. First of 
all, energy may irreversibly be transferred from the system to the environment 
thereby giving rise to the phenomenon of dissipation. In addition, the fluctuating 
force exerted by the environment on the system causes fluctuations of the system 
degree of freedom which manifest itself for example as Brownian motion. While 
these two effects occur both for classical as well as quantum systems, there exists 
a third phenomenon which is specific to the quantum world. As a consequence 
of the entanglement between system and environmental degrees of freedom a 
coherent superposition of quantum states may be destroyed in a process referred 
to as decoherence. This effect is of major concern if one wants to implement a 
quantum computer. Therefore, decoherence is discussed in detail in Chap. 5. 

Quantum computation, however, is by no means the only topic where the 
coupling to an environment is relevant. In fact, virtually no real system can be 
considered as completely isolated from its surroundings. Therefore, the phenom- 
ena listed in the previous paragraph play a role in many areas of physics and 
chemistry and a series of methods has been developed to address this situation. 
Some approaches like the master equations discussed in Chap. 2 are particularly 
well suited if the coupling to the environment is weak, a situation desired in 
quantum computing. On the other hand, in many solid state systems, the envi- 
ronmental coupling can be so strong that weak coupling theories are no longer 
valid. This is the regime where the path integral approach has proven to be very 
useful. 

It would be beyond the scope of this chapter even to attempt to give a com- 
plete overview of the use of path integrals in the description of dissipative quan- 
tum systems. In particular for a two-level system coupled to harmonic oscillator 
degrees of freedom, the so-called spin-boson model, quite a number of approx- 
imations have been developed which are useful in their respective parameter 
regimes. This chapter rather attempts to give an introduction to path integrals 
for readers unfamiliar with but interested in this method and its application to 
dissipative quantum systems. 

In this spirit. Sect. ^ gives an introduction to path integrals. Some aspects 
discussed in this section are not necessarily closely related to the problem of 
dissipative systems. They rather serve to illustrate the path integral approach 
and to convey to the reader the beauty and power of this approach. In Sect. S 
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we elaborate on the general idea of the coupling of a system to an environment. 
The path integral formalism is employed to eliminate the environmental degrees 
of freedom and thus to obtain an effective description of the system degree of 
freedom. The results provide the basis for a discussion of the damped harmonic 
oscillator in Sect. ^. Starting from the partition function we will examine several 
aspects of this dissipative quantum system. 

Readers interested in a more in-depth treatment of the subject of quantum 
dissipation are referred to existing textbooks. In particular, we recommend the 
book by U. Weiss |] which provides an extensive presentation of this topic 
together with a comprehensive list of references. Chapter 4 of may serve as a 
more concise introduction complementary to the present chapter. Path integrals 
are discussed in a whole variety of textbooks with an emphasis either on the 
physical or the mathematical aspects. We only mention the book by H. Kleinert 

which gives a detailed discussion of path integrals and their applications in 
different areas. 



2 Path Integrals 
2.1 Introduction 

The most often used and taught approach to nonrelativistic quantum mechan- 
ics is based on the Schrodinger equation which possesses strong ties with the 
the Hamiltonian formulation of classical mechanics. The nonvanishing Poisson 
brackets between position and momentum in classical mechanics lead us to in- 
troduce noncommuting operators in quantum mechanics. The Hamilton function 
turns into the Hamilton operator, the central object in the Schrodinger equation. 
One of the most important tasks is to find the eigenfunctions of the Hamilton 
operator and the associated eigenvalues. Decomposition of a state into these 
eigenfunctions then allows us to determine its time evolution. 

As an alternative, there exists a formulation of quantum mechanics based 
on the Lagrange formalism of classical mechanics with the action as the cen- 
tral concept. This approach, which was developed by Feynman in the 1940's 
[HP , avoids the use of operators though this does not necessarily mean that the 
solution of quantum mechanical problems becomes simpler. Instead of finding 
eigenfunctions of a Hamiltonian one now has to evaluate a functional integral 
which directly yields the propagator required to determine the dynamics of a 
quantum system. Since the relation between Feynman's formulation and classi- 
cal mechanics is very close, the path integral formalism often has the important 
advantage of providing a much more intuitive approach as we will try to convey 
to the reader in the following sections. 



2.2 Propagator 

In quantum mechanics, one often needs to determine the solution Itpit)) of the 
time-dependent Schrodinger equation 

i^f^ = m^), (1) 
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where H is the Hamiltonian describing the system. Formally, the solution of 
may be written as 



) 1^(0)) . (2) 



Here, the time ordering operator T is required because the operators correspond- 
ing to the Hamiltonian at different times in general due not commute. In the 
following, we will restrict ourselves to time-independent Hamiltonians where 
simplifies to 

\m) ^ exv(^-^H?j\m) ■ (3) 

As the inspection of (||) and (||) demonstrates, the solution of the time-dependent 
Schrodinger equation contains two parts: the initial state |V'(0)) which serves as 
an initial condition and the so-called propagator, an operator which contains all 
information required to determine the time evolution of the system. 
Writing (^) in position representation one finds 

{x\m) = J dx'{x\exp(^-^H?j \x'){x'\i:{0)) (4) 



or 

i^ix, t)= j dx'K{x, t, x', 0)V'(a;', 0) (5) 

with the propagator 

K{x,t,x\0) = {x\exp(^~^H?j \x') . (6) 

It is precisely this propagator which is the central object of Feynman's formu- 
lation of quantum mechanics. Before discussing the path integral representation 
of the propagator, it is therefore useful to take a look at some properties of the 
propagator. 

Instead of performing the time evolution of the state |'0(O)) into in one 

step as was done in equation (^), one could envisage to perform this procedure 
in two steps by first propagating the initial state |^/'(0)) up to an intermediate 
time ti and taking the new state \'ijj{ti)) as initial state for a propagation over 
the time t — ti. This amounts to replacing (|^) by 

\m) = exp(^-^H{t - t,)^ exp(^-^Hh^ |V(0)) (7) 

or equivalently 

*)= J da^y dx"K{x, t, x" , ti)K{x", ti,x', 0)i^{x', 0) . (8) 
Comparing (^) and (|8|), we find the semigroup property of the propagator 

K{x, t, x', 0) = / dx"K{x, t, x", ti)K{x", ti,x\ 0) . (9) 
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Fig. 1. According to tlie semigroup property (^) tlie propagator K{x,t,x' ,0) may be 
decomposed into propagators arriving at some time ti at an intermediate point x" and 
propagators continuing from tliere to tlie final point x 



This result is visualized in Fig. |] where the propagators between space-time 
points are depicted by straight lines connecting the corresponding two points. 
At the intermediate time ti one has to integrate over all positions x" . This insight 
will be of use when we discuss the path integral representation of the propagator 
later on. 

The propagator contains the complete information about the eigenenergies 
En and the corresponding eigenstates \n). Making use of the completeness of the 
eigenstates, one finds from (||) 

if(a;,i,x',0) -^cxpU-i^„tj i^n{x)Mx'y ■ (10) 

Here, the star denotes complex conjugation. Not only does the propagator con- 
tain the eigenenergies and eigenstates, this information may also be extracted 
from it. To this end, we introduce the retarded Green function 

Grix, t, x', 0) = K{x, t, x', 0)&{t) (11) 

where &{t) is the Heavisidc function which equals 1 for positive argument t and 
is zero otherwise. Performing a Fourier transformation, one ends up with the 
spectral representation 

Gr{x,x',E) ^ dtexp(^^Et^ G,{t) 

■lpn{x)tl'n{x^ 



E 



'A* (12) 



E-En+ie ' 



where e is an infinitely small positive quantity. According to (p^), the poles of 
the energy-dependent retarded Green function indicate the eigenenergies while 
the corresponding residua can be factorized into the eigenfunctions at positions 
X and x'. 
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2.3 Free Particle 

An important step towards the path integral formulation of quantum mechanics 
can be made by considering the propagator of a free particle of mass m. The 
eigenstates of the corresponding Hamiltonian 

(13) 

are momentum eigenstates 



with a momentum eigenvalue p out of a continuous spectrum. Inserting these 
eigenstates into the representation (II^) of the propagator, one finds by virtue of 



dxexp(— iax ) = W — = W — expl — i— ] (15) 



for the propagator of the free particle the result 

n2 



K{xut,Xi,Q) = ^ j ^P^^p(^~'^|^^^ expQp(a;f - x-,] 

I m(x[ — Xi)'^ 



2TT\ht \h 2t 



(16) 



It was already noted by Dirac Q that the quantum mechanical propagator 
and the classical properties of a free particle are closely related. In order to 
demonstrate this, we evaluate the action of a particle moving from Xi to Xf in 
time t. From the classical path 

g 

Xclis) = Xi + {xi ~ Xi)- (17) 

obeying the boundary conditions a;ci(0) = Xi and Xci{t) ~ X{, the corresponding 
classical action is found as 

c ™ Z"*^ -2 m{xf-Xi)^ 

Sc^ = yI d.x,, = - (18) 

This result enables us to express the propagator of a free particle entirely in 
terms of the classical action as 

K(xf,i,x,,0)=(^-^ ) exp(^-^ei(xf,i,x,0)j . (19) 

This result is quite remarkable and one might suspect that it is due to a pecu- 
liarity of the free particle. However, since the propagation in a general potential 
(in the absence of delta function contributions) may be decomposed into a series 
of short-time propagations of a free particle, the result ([l^) may indeed be em- 
ployed to construct a representation of the propagator where the classical action 
appears in the exponent. In the prefactor, the action appears in the form shown 



in equation ( |19| ) only within the semiclassical approximation (cf. Sect. 2. 
for potentials where this approximation turns out to be exact. 



6 Gert-Ludwig Ingold 



2.4 Path Integral Representation of Quantum Mechanics 

While avoiding to go too deeply into the mathematical details, we nevertheless 
want to sketch the derivation of the path integral representation of the propaga- 
tor. The main idea is to decompose the time evolution over a finite time t into 
slices of short time intervals At = t/N where we will eventually take the limit 
N — > oo. Denoting the operator of the kinetic and potential energy by T and V, 
respectively, we thus find 



exp(^-^HtJ ^ exp(^-^{T + V)A?j 



-I N 

(20) 



For simplicity, we will assume that the Hamiltonian is time-independent even 
though the following derivation may be generalized to the time-dependent case. 
We now would like to decompose the short-time propagator in (|2(]| ) into a part 
depending on the kinetic energy and another part containing the potential en- 
ergy. However, since the two operators do not commute, we have to exercise 
some caution. From an expansion of the Baker- Hausdorfi' formula one finds 



exp(^-^{T + V)A?j ^exp(^-^TA?j exp(^-^VA?j + ^[T,V]{Atf 



(21) 



where terms of order {At)^ and higher have been neglected. Since we are inter- 
ested in the limit At — > 0, we may neglect the contribution of the commutator 
and arrive at the Trotter formula 

exp(^-^(r + V)t^ = \im^ [UiAt)f (22) 

with the short time evolution operator 

U{At) = cxp(^- ^TAt^ exp(^-^VAt^ . (23) 

What we have presented here is, of course, at best a motivation and certainly 
does not constitute a mathematical proof. We refer readers interested in the 
details of the proof and the conditions under which the Trotter formula holds to 
the literature |7|. 

In position representation one now obtains for the propagator 



K{xf, t, Xi, 0) = lim 



l-oo (n-1 \ 

j ( n dx, 1 (xf |C/(Z\t)|a;Ar_ 



(24) 

X {xi\U{At)\xi) . 



Since the potential is diagonal in position representation, one obtains together 
with the expression (|l|) for the propagator of the free particle for the matrix 
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element 



2-K\hAt 



exp 



m {xj+i — Xj)'^ 



exp[ ~-V{xj)At 
-V{xj)At 



At 



(25) 



We thus arrive at our final version of the propagator 



K{x[, t, Xi, 0) — lim 




(26) 



X exp 



where xq and xm should be identified with Xi and Xf , respectively. The discretiza- 
tion of the propagator used in this expression is a consequence of the form ( pl[ ) 
of the Baker-Hausdorff relation. In lowest order in we could have used a 
different decomposition which would have led to a different discretization of the 
propagator. For a discussion of the mathematical subtleties we refer the reader 
to@. 

Remarking that the exponent in ( |2^ ) contains a discretized version of the 
action 



S[x]= As{^x^-V{x)) 



we can write this result in short notation as 



K{xi,t,Xi,Q) 



I?xexp( 



(27) 



(28) 



The action ( pTj ) is a functional which takes as argument a function x[s) and 
returns a number, the action S[x\. The integral in ( p8[ ) therefore is a functional 
integral where one has to integrate over all functions satisfying the boundary 
conditions a:(0) = Xi and x{t) — Xf. Since these functions represent paths, one 
refers to this kind of functional integrals also as path integral. 

The three lines shown in Fig. y represent the infinity of paths satisfying 
the boundary conditions. Among them the thicker line indicates a special path 
corresponding to an extremum of the action. According to the principal of least 
action such a path is a solution of the classical equation of motion. It should be 
noted, however, that even though sometimes there exists a unique extremum, in 
general there may be more than one or even none. A demonstration of this fact 
will be provided in Sect. 2.7 where we will discuss the driven harmonic oscillator. 

The other paths depicted in Fig. |^ may be interpreted as quantum fluctu- 
ations around the classical path. As we will see in Sect. 2.S, the amplitude of 
these fluctuations is typically of the order of ^/h. In the classical limit ft ^ 
therefore only the classical paths survive as one should expect. 
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t 

Fig. 2. The thick hne represents a classical path satisfying the boundary conditions. 
The thinner lines are no solutions of the classical equation of motion and may be 
associated with quantum fluctuations 

Before explicitly evaluating a path integral, we want to discuss two examples 
which will give us some insight into the difference of the approaches offered by 
the Schrodinger and Feynman formulation of quantum mechanics. 



2.5 Particle on a Ring 

We confine a particle of mass rn to a ring of radius R and denote its angular 
degree of freedom by (j). This system is described by the Hamiltonian 

"^'2^W ^^^^ 

Requiring the wave function to be continuous and differentiable, one finds the 
stationary states 

(</>) = ^exp(i^(/)) (30) 
with £ = 0,±1,±2,... and the eigenenergies 

These solutions of the time-independent Schrodinger equation allow us to con- 
struct the propagator 

K{<PutAuO) = ^ E exp(i£(0f-0O-i^;;^O • (32) 

We now want to derive this result within the path integral formalism. To this 
end we will employ the propagator of the free particle. However, an important 
difference between a free particle and a particle on a ring deserves our attention. 
Due to the ring topology we have to identify all angles (j) + 27rn, where n is an 
integer, with the angle (f). As a consequence, there exist infinitely many classical 
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n — n — 1 

Fig. 3. On a ring, the angles (pf and (pi + 27rn have to be identified. As a consequence, 
there exist infinitely many classical paths connecting two points on a ring, which may 
be identified by their winding number n 

paths connecting (j>i and (pi. All these paths are topologically different and can 
be characterized by their winding number n. As an example. Fig. ^ shows a path 
for n = and n = 1. Due to their different topology, these two paths (and any 
two paths corresponding to different winding numbers) cannot be continuously 
transformed into each other. This implies that adding a fluctuation to one of the 
classical paths will never change its winding number. 

Therefore, we have to sum over all winding numbers in order to account for 
all possible paths. The propagator thus consists of a sum over free propagators 
corresponding to different winding numbers 

i^fA A n\ D / /i TOi?2 (0f _ 0; _ 27rn)2\ 

K(0M,0i,o)^ E ^V^'^^ft^ 1 j • 

n— — oo ^ ' 

Here, the factor R accounts for the fact that, in contrast to the free particle, the 
coordinate is given by an angle instead of a position. 

The propagator ( p3| ) is 27r-periodic in 0f — 0; and can therefore be expressed 
in terms of a Fourier series 

oo 

K{(j)i, 0i, 0) = C£ exp [ie{(j)i - 0i)] . (34) 

i=-oo 

The Fourier coefficients are found to read 

1 / \ 
c^ = ^exp(^-:^tj 

which proves the equivalence of ( |3^ ) with our previous result (|3^). We thus 
have obtained the propagator of a free particle on a ring both by solving the 
Schrodinger equation and by employing path integral methods. These two ap- 
proaches make use of complementary representations. In the first case, this is 
the angular momentum representation while in the second case, one works in 
the phase representation and sums over winding numbers. 




(35) 
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Fig. 4. The reflection at the waUs of a box leads to an infinite number of possible 
trajectories connecting two points in the box 



2.6 Particle in a Box 

Another textbook example in standard quantum mechanics is the particle in a 
box of length L confined by infinitely high walls at a; = and x = L. From the 
eigenvalues 

with j — 1,2,... and the corresponding eigenfunctions 



(a;) = Xlj^ sin (^TTj - j (37) 



the propagator is immediately obtained as 



if(xf,i,Xi,0) = - ^expf-i^^-^n sin(^7rj^^ sin(^7rj^) . (38) 



It took some time until this problem was solved within the path integral 
approach |H,|lO). Here, we have to consider all paths connecting the points Xi and 
Xi within a period of time t. Due to the reflecting walls, there again exist infinitely 
many classical paths, five of which are depicted in Fig. ^. However, in contrast to 
the case of a particle on a ring, these paths are no longer topologically distinct. 
As a consequence, we may deform a classical path continuously to obtain one of 
the other classical paths. 

If, for the moment, we disregard the details of the reflections at the wall, the 
motion of the particle in a box is equivalent to the motion of a free particle. The 
fact that paths are folded back on themselves can be accounted for by taking 
into account replicas of the box as shown in Fig. |^. Now, the path does not 
necessarily end at xf' — X{ but at one of the mirror images where n is an 
arbitrary integer. In order to obtain the propagator, we will have to sum over 
all replicas. Due to the different geometry we need to distinguish between those 
paths arising from an even and an odd number of reflections. From Fig. ^ one 
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Fig. 5. Instead of a particle getting reflected at the walls of the box one may think of 
a free particle moving from the starting point in the box to the end point in one of the 
replicas of the box 

can see that for an odd number 2n— 1 of reflections, the end point Ues at 2nL — X{ 
and the contribution to the full propagator therefore is given by 

On the other hand, for an even number 2n of reflections, the end point is located 
at 2nL + X{ and wc find 

io„\ , , I m / i 'm(2nL + Xf — x;)'^\ 

- i2^tHli- 2^^^) ■ (40) 

However, it is not obvious that just summing up the propagators ( |39| ) and ( pO[ ) 
for all n will do the job. 

In order to clarify this point, we start with the somewhat simpler situation 
of just one wall and take a look at all paths running between xi and Xf in time 
t. As can be seen from the space-time diagram in Fig. ^ there are paths which 
do not cross the wall and which therefore contribute to the path integral. On 
the other hand, there exist also paths which cross the wall an even number of 
times. Since these paths spend some time in the forbidden region, they do not 
contribute to the path integral. 

It requires some thinking to ensure that only paths not crossing the wall are 
taken into account. Our strategy will consist in first writing down a propagator 
Kfy-ee which disregards the wall. Then, we have to subtract off the contributions 
of all the paths which cross the wall. This can be done by constructing a path 
with the same action as the original path. To this end we take the original path 
up to the last crossing with the wall and then continue along the mirror image 
of the original path. We thus end up at the mirror image —Xf of the original end 
point xi. Note that a path running from Xi to —X{ necessarily crosses the wall 
at least once. As a consequence, subtracting the propagator between these two 
points eliminates all original paths which do not remain in the region x > 0. 
We therefore obtain our desired result, the propagator ifwaii in the presence of 
a wall, by subtracting a propagator going to the refiected end point from the 
unconstrained propagator to the orig inal end point 



K„an{xi, t, Xi, 0) = Kfrcc{Xf, t, Xi, 0) - i^trco(-a;f , t, Xi, 0) 



(41) 
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Fig. 6. A path crossing the wall is cancelled by a path running to the mirror point of 
the end point 

This result bears much resemblance with the method of image charges in elec- 
trostatics. After giving it some thought, this should not be too surprising since 
the free Schrodinger equation and the Poisson equation are formally equivalent. 
According to the method of image charges one may account for a metallic plate 
(i.e. the wall) by putting a negative charge (i.e. the mirrored end point) com- 
plementing the positive charge (i.e. the original end point). For the propagator 
this results in the difference appearing in (^). 

Let us now come back to our infinitely deep potential well with two walls. This 
problem corresponds to the electrostatics of a charge between two parallel metal 
plates. In this case, the method of image charges leads to an infinite number 
of charges of alternating signs. The original positive charge gives rise to two 
negative charges which are each an image corresponding to one of the two metal 
plates. In addition, however, these images have mirror images corresponding to 
the other metal plate and this process has to be carried on ad infinitum. 

Expressing the propagator of the particle in the box in terms of the free 
propagator works in exactly the same way. A path intersecting both walls is 
subtracted twice, i.e. one time too often. Therefore, one contribution has to be 
restored which is done by adding another end point. Continuing the procedure 
one ends up with an infinite number of end points, some of which we have shown 
m Fig. |. As a consequence, we can attribute a sign to each end point in this 
figure. The general rule which follows from these considerations is that each 
reflection at a wall leads to factor —1. The propagator therefore can be written 
as 




(42) 
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The symmetries 

K{x[ + 2L, t, Xi, 0) = K{xf, t, Xi, 0) 
K{—Xi, t, Xi, 0) = ~K{xi, t, Xi, 0) 

suggest to expand the propagator into the Fourier series 



K{xut,X{,Qi) = ^aj(a;i,t) sin (^^rj-^ 



Its Fourier coefRcients are obtained from (42) as 

r-L 



aj(xi,t) = -^ / dxifim(^TTj^^ K{xi,t,xuQ) 

J — Z/ 



2 

^-sml 



(43) 
(44) 



(45) 



(46) 



where the energies Ej are the eigenenergies of the box defined in (js^) • Inserting 
(Eel) into (EBh we thus recover our previous resuh (^. 



2.7 Driven Harmonic Oscillator 

Even though the situations deah with in the previous two sections have been con- 
ceptually quite interesting, we could in both cases avoid the explicit calculation 
of a path integral. In the present section, we will introduce the basic techniques 
needed to evaluate path integrals. As an example, we will consider the driven 
harmonic oscillator which is simple enough to allow for an exact solution. In 
addition, the propagator will be of use in the discussion of damped quantum 
systems in later sections. 

Our starting point is the Lagrangian 

i=|i2„™^2^2^^^(,) (47) 

of a harmonic oscillator with mass m and frequency lo. The force /(t) may be 
due to an external field, e.g. an electric field coupling via dipole interaction to a 
charged particle. In the context of dissipative quantum mechanics, the harmonic 
oscillator could represent a degree of freedom of the environment under the 
influence of a force exerted by the system. 

According to ( p8| ) we obtain the propagator K{xi^ t, xj, 0) by calculating the 
action for all possible paths starting at time zero at x; and ending at time t at 
Xi- It is convenient to decompose the paths 

a;(s) =a;ci(s)+C(s) (48) 

into the classical path Xc\ satisfying the boundary conditions a;ci(0) = Xi, Xci(i) = 
Xi and a fluctuating part ^ vanishing at the boundaries, i.e. ^(0) = ^(<) = 0. The 
classical path has to satisfy the equation of motion 



mXci + vnijp'XcX — f{s) 



(49) 
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obtained from the Lagrangian (p7|). 

For an exactly solvable problem like the driven harmonic oscillator, we could 
replace Xd by any path satisfying x{0) = Xi, x{t) = X{. We leave it as an 
exercise to the reader to perform the following calculation with Xci{s) of the 
driven harmonic oscillator replaced by xi + {x{ — Xi)s/t. However, it is important 



to note that within the semiclassical approximation discussed in Sect. 2. 
expansion around the classical path is essential since this path leads to the 
dominant contribution to the path integral. 
With (148) we obtain for the action 



S = j\s{^x^~^oj^x^+xf{s)] 

( Y'^ci - + a;ci/(s)) + ds (mxdi - muj'^XciC + ^f{s)j 



(50) 

For our case of a harmonic potential, the third term is independent of the bound- 
ary values Xi and X{ as well as of the external driving. The second term vanishes 
as a consequence of the expansion around the classical path. This can be seen 
by partial integration and by making use of the fact that Xd is a solution of the 
classical equation of motion: 

ds(^mXciS,—mijj'^XciS,+£,f{s))=— / ds(mXci + muj'^Xci — f{s))£, — . (51) 

We now proceed in two steps by first determining the contribution of the 
classical path and then addressing the fluctuations. The solution of the classical 
equation of motion satisfying the boundary conditions reads 

, , sinfws) sin(cj(t — s)) , 
Xciis) = Xf + Xi \ ) 52) 
sm(wr) sm(wi) 



1 



dusm{uj(s — u))f(u) ^-4 — r / dusm{uj(t ~ u))f(u) 

sm{Lut) 7o 



A peculiarity of the harmonic oscillator in the absence of driving is the ap- 
pearance of conjugate points at times Tn — (7r/a;)n where n is an arbitrary 
integer. Since the frequency of the oscillations is independent of the amplitude, 
the position of the oscillator at these times is determined by the initial position: 
x{T2n+i) = — a^i and x(T2„) = Xj. This also illustrates the fact mentioned on 
p. 0, that depending on the boundary conditions there may be more than one 
or no classical solution. 
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Fig. 7. In a harmonic potential all trajectories emerging from the same starting point 
converge at conjugate points at multiples of half an oscillation period 



The task of evaluating the action of the classical path may be simplified by 
a partial integration 



/ 771/ Tfl \ 

- y ds y—Xc\Xc\ + l^^'^^lx - Xcif{s)j (53) 





m 



Tn 1 f 

^ {x{Xci{t) - XiXci{0)) + - J dsxciis)f{s) 

where we have made use of the classical equation of motion to obtain the third 
line. From the solution (|5^) of the classical equation of motion we get 

, , Xf — XiCOS(ujt) 1 , , ^^ ^ . ^ 

x^^(0) = uj- — . ' ; ' dssm{uj{t- s))f{s) (54) 



sin{Lut) m sin{ujt) 



Xc\{t) = LO 



Xi cos(ujt) — Xi 



1 



ds sin(a;s)/(s) . 



sa\{uot) m s\x\{uit) Jq 

Inserting initial and final velocity into (p3|) we find for the classical action 



(55) 



2 sin(w<) 

Xl 



[{x? + Xf) cos(cji) — 2xiXf] 



+ — — - / dssin(ws)/(s) + — A— / ds sin(a;(t - s))/(s) (56) 
sm(wt) 7o sm{ujt) J„ ^ ^ 



1 



ds / du sin(a;it) sin(cj(t — s))/(s)/(w) 



ma;sin(a;t) Jq Jq 

As a second step, we have to evaluate the contribution of the fluctuations 
which is determined by the third term in (|5^). After partial integration this term 
becomes 

(57) 











ds2 



16 Gert-Ludwig Ingold 



Here, the superscript '(2)' indicates that this term corresponds to the contri- 
bution of second order in ^. In view of the right-hand side it is appropriate to 
expand the fluctuation 



= ^a„e«(5) (58) 



into eigenfunctions of 



d^ 
ds2 



e« = A„e« (59) 



with ^n(O) = S,n{t) — 0. As eigenfunctions of a selfadjoint operator, the ^„ are 
complete and may be chosen orthonormal. Solving (|5^) yields the eigenfunctions 

Cn(.s) = yfsin(7rn^) (60) 

and corresponding eigenvalues 

A„ = -(^)'+t.2. (61) 

We emphasize that (^^ is not the usual Fourier series on an interval of length 
t. Such an expansion could be used in the form 

- [an (cos(27rn^) - l) + &„ sin(27rn^)] (62) 

n— 1 

which ensures that the fluctuations vanish at the boundaries. We invite the 
reader to redo the following calculation with the expansion (^) replacing (|58|). 
While at the end the same propagator should be found, it will become clear why 
the expansion in terms of eigenfunctions satisfying (|5^) is preferable. 

The integration over the fluctuations now becomes an integration over the 
expansion coefficients a„. Inserting the expansion (Kq) into the action one finds 



n=l n=l ^ ^ 

As this result shows, the classical action is only an extremum of the action but 
not necessarily a minimum although this is the case for short time intervals 
t < tt/lo. The existence of conjugate points at times r„ = mr/uj mentioned 
above manifests itself here as vanishing of the eigenvalue A„ . Then the action is 
independent of a„ which implies that for a time interval r„ all paths x^-i + flnCn 
with arbitrary coefficient a„ are solutions of the classical equation of motion. 

After expansion of the fiuctuations in terms of the eigenfunctions (|6^), the 
propagator takes the form 
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In principle, we need to know the Jacobi determinant of the transformation from 
the path integral to the integral over the Fourier coefficients. However, since this 
Jacobi determinant is independent of the oscillator frequency w, we may also 
compare with the free particle. Evaluating the Gaussian fluctuation integrals, 
we find for the ratio between the prefactors of the propagators K^^ and Kq of 
the harmonic oscillator and the free particle, respectively, 

KoeM-ii/h)S,,,o] D ■ ^ ' 

Here, we have introduced the fluctuation determinants for the harmonic oscilla- 
tor 

D^det( — +^'\=l[Xr. (66) 

^ ^ n— 1 

and the free particle 

^°=^^Kd^) = n<- (67) 

The eigenvalues for the free particle 



n=l 



are obtained from the eigenvalues (|6^) of the harmonic oscillator simply by 
setting the frequency uj equal to zero. With the prefactor of the propagator of 
the free particle 

and (|65|), the propagator of the harmonic oscillator becomes 

Kixut,x,,0) = .r^Me^pfU^ . (70) 



For readers unfamiliar with the concept of determinants of differential opera- 
tors we mention that we may define matrix elements of an operator by projection 
onto a basis as is familiar from standard quantum mechanics. The operator rep- 
resented in its eigenbasis yields a diagonal matrix with the eigenvalues on the 
diagonal. Then, as for finite dimensional matrices, the determinant is the product 
of these eigenvalues. 



Each of the determinants (pq) and (67) by itself diver ges. However, we are 



interested in the ratio between them which is well-defined |]12| 

(71) 
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Inserting this result into (^0|) leads to the propagator of the driven harmonic 
oscillator in its final form 



with the classical action defined in (p6[). The Morse index n in the phase factor is 
given by the integer part of Lot/?:. This phase accounts for the changes in sign of 
the sine function |p3| . Here, one might argue that it is not obvious which sign of 
the square root one has to take. However, the semigroup property allows to 
construct propagators across conjugate points by joining propagators for shorter 
time intervals. In this way, the sign may be determined unambiguously |l4| . 

It is interesting to note that the phase factor exp(— in7r/2) in (|7^ ) implies that 
K(x{,2TT/uj,Xi,0) = —K{x[,0,Xi,0) = —S{x[ — Xi), i.e. the wave function after 
one period of oscillation differs from the original wave function by a factor —1. 
The oscillator thus returns to its original state only after two periods very much 
like a spin-1/2 particle which picks up a sign under rotation by 2t: and returns 
to its original state only after a 47r-rotation. This effect might be observed in 
the case of the harmonic oscillator by letting interfere the wave functions of two 
oscillators with different frequency [T^ . 

2.8 Semiclassical Approximation 

The systems considered so far have been special in the sense that an exact 
expression for the propagator could be obtained. This is a consequence of the 
fact that the potential was at most quadratic in the coordinate. Unfortunately, 
in most cases of interest the potential is more complicated and apart from a few 
exceptions an exact evaluation of the path integral turns out to be impossible. 
To cope with such situations, approximation schemes have been devised. In the 
following, we will restrict ourselves to the most important approximation which 
is valid whenever the quantum fluctuations are small or, equivalently, when the 
actions involved are large compared to Planck's constant so that the latter may 
be considered to be small. 

The decomposition of a general path into the classical path and fluctuations 
around it as employed in (^) in the previous section was merely a matter of 
convenience. For the exactly solvable case of a driven harmonic oscillator it 
is not really relevant how we express a general path satisfying the boundary 
conditions. Within the semiclassical approximation, however, it is decisive to 
expand around the path leading to the dominant contribution, i.e. the classical 
path. From a more mathematical point of view, we have to evaluate a path 
integral over exp(iiS'/?i) for small H. This can be done in a systematic way by the 
method of stationary phase where the exponent has to be expanded around the 
extrema of the action S. 




(72) 
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Fig. 8. In stationary phase approximation only a small region around the extremum 
contributes to the integral. For the example shown here, the extremum lies at a; = 



At this point it may be useful to give a brief reminder of the method of 
stationary phase. Suppose we want to evaluate the integral 



1(a) 



dxg{x) exp(ia/(x)) 



(73) 



in the limit of very large a. Inspection of Fig. ||, where f{x) = x^, suggests that 
the dominant contribution to the integral comes from a region, in our example 
of size 1/v^j around the extremal (or stationary) point of the function f{x). 
Outside of this region, the integrand is rapidly oscillating and therefore gives to 
leading order a negligible contribution. Since for large a, the region determining 
the integral is very small, we may expand the function /(x) locally around the 
extremum xq 

fix) « f{xo) + lf"{xo){x - xo)^ + ... (74) 

and replace g{x) by g{xo). Neglecting higher order terms, which is allowed if 
f"{xo) is of order one, we are left with the Gaussian integral 



(75) 



/(a) « g(a;o)exp(ia/(xo)) / da;exp( -/"(xo)(a; - xq)' 



27r r TT 

|/"(xo)| ^^^°^ [i"/(2^o) + i-sgn(/"(a;o)) 



where sgn(/"(xo)) denotes the sign of /"(xq). If f{x) possesses more than one 
extremum, one has to sum over the contributions of all extrema unless one 
extremum can be shown to be dominant. 

We now apply the stationary phase approximation to path integrals where 
l/h plays the role of the large parameter. Since the action is stationary at clas- 
sical paths, we are obliged to express the general path as 



x{s) = a;ci(s) + $(s) 



(76) 



where Xd is the classical path (or one of several possible paths) satisfying the 
boundary conditions and ^ represents the fluctuations around the classical path. 
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With this decomposition the action becomes 



S 



ds (yici - V{x,{)) + y ds (mxci " V'{x,i)i) (77) 
t 



As[^s,'-\v"{x.{)e] + 



It is instructive to compare this result with the action (|5C|) for the driven har- 
monic osciUator. Again, the first term represents the classical action. The second 
term vanishes as was shown explicitly in for the driven oscillator. In the 
general case, one can convince oneself by partial integration of the kinetic part 
and comparison with the classical equation of motion that this term vanishes 
again. This is of course a consequence of the fact that the classical path, around 
which we expand, corresponds to an extremum of the action. The third term on 
the right-hand-side of (|7^) is the leading order term in the fluctuations as was 
the case in (|50|). There is however an important difference since for anharmonic 
potentials the second derivative of the potential V" is not constant and there- 
fore the contribution of the fluctuations depends on the classical path. Finally, 
in general there will be higher order terms in the fluctuations as indicated by the 
dots in (fr^ ) . The semiclassical approximation consists in neglecting these higher 
order terms so that after a partial integration, we get for the action 



(78) 



where the index 'sc' indicates the semiclassical approximation. 

Before deriving the propagator in semiclassical approximation, we have to 
discuss the regime of validity of this approximation. Since the first term in ( [78| ) 
gives only rise to a global phase factor, it is the second term which determines the 
magnitude of the quantum fluctuations. For this term to contribute, we should 
have /h ^ 1 so that the magnitude of typical fluctuations is at most of order 
^/h. The term of third order in the fluctuations is already smaller than the second 
order term by a factor {^/h)'^/h = \Jli. If Planck's constant can be considered 
to be small, we may indeed neglect the fluctuation contributions of higher than 
second order except for one exception: It may happen that the second order term 
does not contribute, as has been the case at the conjugate points for the driven 



harmonic oscillator in Sect. 2/7. Then, the leading nonvanishing contribution 
becomes dominant. For the following discussion, we will not consider this latter 
case. 



In analogy to Sect. 2.7 we obtain for the propagator in semiclassical approx- 
imation 

K{x,X^u^) = J^a\^^^v{\s,^ (79) 
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where 



(80) 



and Dq is the fluctuation determinant ( pT]) of the free particle. 

Even though it may seem that determining the prefactor is a formidable task 
since the fluctuation determinant for a given potential has to be evaluated, this 
task can be greatly simplified. In addition, the following considerations offer the 
benefit of providing a physical interpretation of the prefactor. In our evaluation 
of the prefactor we follow Marinov The main idea is to make use of the 
semigroup property (^) of the propagator 



C{x{, t, Xi, 0) exp 



:S'ci(xf,t, XijO) 



(81) 



= / dx C[xi,t,x ,t')C{x' ,t' ,Xi,Q) eyi\) 



[S'ci(a;f , t, X , t') + S'ci(x', t' , Xi, 0)] 



where the prefactor C depends on the fluctuation contribution. We now have 
to evaluate the x'-integral within the semiclassical approximation. According to 
the stationary phase requirement discussed above, the dominant contribution to 
the integral comes from x' = a;o(a;f , Xj, i, t') satisfying 



dS,i{xut,x\t') 



dx' 



95ci(x',t',Xi,0) 



X' —Xq 



dx' 



. 



(82) 



X' —Xq 



According to classical mechanics these derivatives are related to initial and final 
momentum by WJl 



dxi 



-Pi 



Xf,t{ ,ti 



dxi 



Pi 



SO that (82) can expressed as 



p{t' - e) ^ p(t' + e) 



(83) 



(84) 



The point xq thus has to be chosen such that the two partial classical paths 
can be joined with a continuous momentum. Together they therefore yield the 
complete classical path and in particular 



5'ci(a;f , t, Xi, 0) = Sc\(xu t, xq, t') + S'ci(a:o, t' , Xi, 0) 



(85) 



This relation ensures that the phase factors depending on the classical actions 
on both sides of ( |8l| ) are equal. 

After having identified the stationary path, we have to evaluate the integral 
over x' in (^). Within semiclassical approximation this Gaussian integral leads 
to 

C{xf,t,xuO) 



C{xi, t, xo, t')C{xo, t', Xi, 0) 



2TTih Oxq 



Q1 X-l/2 
[S'cl(Xf , t, Xo,t') + Scl{xo, t', Xi, 0)] 
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In order to make progress, it is useful to take the derivative of ([S^) with respect 
to X{ and xi. Keeping in mind that xq depends on these two variables one finds 



d'^Sc\{xi,t,Xi,0) _ d'^Sci{x{,t,xojt') dxo d'^Sciixo,t' ,Xi,0) dxo 



dxfdxi 



dxfdxo dxi 



dxidxQ dx[ 



+ -Q^[^ci{x{,t,Xo,t ) + i)ci(xo,i '2:i,t)J — — 



(87) 



dxi dx{ 



Similarly, one finds by taking derivatives of the stationary phase condition ( p2[ ) 

— Sci{x{,t,xo,t') 

(88) 



dxp 
dx{ 



dxfXo 



92 



[S'cl(Xf , t, Xo, t') + Sclixo,t',Xi, 0)] 



and 



dxi 



dxiXQ 



S'ci(xo,t',a;i,0) 



[Scl{X[,t, Xq, t') + Scl{xo,t', Xi, 0)] 



(89) 



dxi 



These expressions allow to eliminate the partial derivatives of xq with respect 
to Xi and X{ appearing in (|8^) and one finally obtains 



■ [Scl{x[,t, Xq, t') + S'cl(xo, t', Xi, 0)] 



(90) 



92 



dxidxf 



S'ci(xf,t, a;i,0) 



d'^Sci{xf,t,xo,t') d'^Sc\{xo,t',Xi,0) 



dxfdxo 



dxidxo 



Inserting this result into (p6[), the prefactor can be identified as the so-called 
Van Vleck-Pauli-Morette determinant lT§,|ir 



C{xf,t,Xi,0) 



I 



2nih 



d^Sciixf,t,Xi,0) 
dxfdxi 



1/2 



(91) 



so that the propagator in semiclassical approximation finally reads 
K{xf,t,xi,0) (92) 



1 

2Trh 



d'^Sciixf,t,Xi,0) 
dxfdxi 



1/2 



exp 



where the Morse index n denotes the number of sign changes of Sd/ dxfdxi 
[ p!3| . We had encountered such a phase factor before in the propagator ( [7^ of 
the harmonic oscillator. 
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As we have already mentioned above, derivatives of the action with respect 
to position are related to momenta. This allows to give a physical interpretation 
of the prefactor of the propagator as the change of the end point of the path as 
a function of the initial momentum 



dxidxf J dpi 

A zero of this expression, or equivalently a divergence of the prefactor of the 
propagator, indicates a conjugate point where the end point does not depend on 
the initial momentum. 

To close this section, we want to compare the semiclassical result (^2|) with 
exact results for the free particle and the harmonic oscillator. In our discussion 



of the free particle in Sect. 2.3 we already mentioned that the propagator can 
be expressed entirely in terms of classical quantities. Indeed, the expression ( pj[ ) 
for the propagator of the free particle agrees with (|9^). 

For the harmonic oscillator, we know from Sect. ^.7| that the prefactor does 
not depend on a possibly present external force. We may therefore consider the 



action (56) in the absence of driving /(s) = which then reads 



= o ■ / ^^ [i^f + ) cos(wi) - 2xiXi] . (94) 
Taking the derivative with respect to xi and Xf one finds for the prefactor 

/ \ 1/2 

Cixut,xuO)=i ^ . , ,J 95 



which is identical with the prefactor in our previous result ([72|). As expected, 
for potentials at most quadratic in the coordinate the semiclassical propagator 
agrees with the exact expression. 



2.9 Imaginary Time Path Integral 

In the discussion of dissipative systems we will be dealing with a system cou- 
pled to a large number of environmental degrees of freedom. In most cases, the 
environment will act like a large heat bath characterized by a temperature T. 
The state of the environment will therefore be given by an equilibrium density 
matrix. Occasionally, we may also be interested in the equilibrium density ma- 
trix of the system itself. Such a state may be reached after equilibration due to 
weak coupling with a heat bath. 

In order to describe such thermal equilibrium states and the dynamics of 
the system on a unique footing, it is desirable to express equilibrium density 
matrices in terms of path integrals. This is indeed possible as one recognizes by 
writing the equilibrium density operator in position representation 



Pfs{x,x') = — (a;|exp(-/3i?)|x') 



(96) 
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with the partition function 



da;(x| exp(— /3i/)|a;) . 
Comparing with the propagator in position representation 
K{x,t,x',0) = {x\e^p(-^Ht] \x') 



(97) 



(98) 



one concludes that apart from the partition function the equihbrium density 
matrix is equivalent to a propagator in imaginary time t = —ih[3. 

After the substitution a = is the action in imaginary time —ihp reads 



(99) 



Here and in the following, we use greek letters to indicate imaginary times. 
Motivated by the right-hand side of (p9[) we define the so-called Euclidean action 





m 


(Ax\ 


/ ds 






/o 









m 




= i / dcr 




m 


Jq 


y 





S^[x] 



up 



dcr 



-x^ + V{x) 



(100) 



Even though one might fear a lack of intuition for motion in imaginary time, this 
results shows that it can simply be thought of as motion in the inverted potential 
in real time. With the Euclidean action ( |10C| ) we now obtain as an important 
result the path integral expression for the (unnormalized) equilibrium density 
matrix 



(x|exp(-/377)|x') 



rx(hfj)=x / -■ \ 

Pxexp -tS^\A ■ (101) 

'x(0)=x' V / 

This kind of functional integral was discussed as early as 1923 by Wiener |2l) in 
the context of classical Brownian motion. 

As an example we consider the (undriven) harmonic oscillator. There is ac- 
tually no need to evaluate a path integral since we know already from Sect. 2.7 
the propagator 



K[xi,t,xi,^) 



27rifi,sin(cjt) 



exp 



. muj {xf + ) cos(ti)t) — 2xiXf 



2h 



sm(Lut) 



(102) 



Transforming the propagator into imaginary time t —> —ihp and renaming Xi 
and X{ into x' and x, respectively, one obtains the equilibrium density matrix 



pp{x,x') 



(103) 



■ exp 



muj {x^ + x'"^) cosh(fi,/3w) — 2xx' 



2h 



smh{hf3uj) 



Z\l 27rhsmh{hpLu) 
The partition function is obtained by performing the trace as 

Z = I dx{x\eM-m\^) = J—^^ 



ih(ft/3w/2) 



(104) 
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which agrees with the expression 

oo 

Z = ^ exp 



-(3huj n 



(105) 



based on the energy levels of the harmonic oscillator. 

Since the partition function often serves as a starting point for the calculation 
of thermodynamic properties, it is instructive to take a closer at how this quan- 
tity may be obtained within the path integral formalism. A possible approach 
is the one we just have sketched. By means of an imaginary time path integral 
one first calculates (x| exp{—f3H)\x) which is proportional to the probability to 
find the system at position x. Subsequent integration over coordinate space then 
yields the partition function. 

However, the partition function may also be determined in one step. To this 
end, we expand around the periodic trajectory with extremal Euclidean action 
which in our case is given by x{a) = 0. Any deviation will increase both the 
kinetic and potential energy and thus increase the Euclidean action. All other 
trajectories contributing to the partition function are generated by a Fourier 
series on the imaginary time interval from to hp 



X{(7) = 



1 



flo 



V2E 



'an cos(i^„cr) + 5„ sin(zy„cr)) 



where we have introduced the so-called Matsubara frequencies 

27r 



(106) 



(107) 



This ansatz should be compared with (|6^) for the fluctuations where oq was 
fixed because the fluctuations had to vanish at the boundaries. For the partition 
function this requirement is dropped since we have to integrate o ver a ll periodic 
trajectories. Furthermore, we note that indeed with the ansatz (106) only the 
periodic trajectories contribute. All other paths cost an infinite amount of action 
due to the jump at the boundary as we will see shortly. 

Inserting the Fourier expansion (10£) into the Euclidean action of the har- 
monic oscillator 



we find 



m 



d(T— (i^ +uj'^x'^) 



(108) 



(109) 



As in Sect. 2.7 we do not want to go into the mathematical details of integration 
measures and Jacobi determinants. Unfortunately, the free particle cannot serve 
as a reference here because its partition function does not exist. We therefore 
content ourselves with remarking that because of 



1 



hf3 



1 



E„.ii^^2sinh(n/3c./2) 



(110) 
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the result of the Gaussian integral over the Fourier coefficients yields the parti- 
tion function up to a frequency independent factor. This enables us to determine 
the partition function in more complicated cases by proceeding as above and us- 
ing the partition function harmonic oscillator as a reference. 

Returning to the density matrix of the harmonic oscillator we finally obtain 



by inserting the partition function (104) into the expression (103) for the density 
matrix 



Pl3(,X,x') =4 



' muj cosh{h(3uj) — 1 



X exp 



smh{h(3uj) 

muj (x^ 



(111) 



x'"^) cosh{h(3u!) — 2xx' 



2h 



smh{h(3uj) 



Without path integrals, this result would require the evaluation of sums over 
Hermite polynomials. 



The expression for the density matrix (111) can be verified in the limits of 
high and zero temperature. In the classical limit of very high temperatures, the 
probability distribution in real space is given by 



2tt 



exp —(3 



exp[-f3V{x) 



(112) 



We thus have obtained the Boltzmann distribution which depends only on the 
potential energy. The fact that the kinetic energy does not play a role can easily 
be understood in terms of the path integral formalism. Excursions in a very short 
time hp cost too much action and are therefore strongly suppressed. 

In the opposite limit of zero temperature the density matrix factorizes into 
a product of ground state wave functions of the harmonic oscillator 



lim Pf3{x, x') ~ 



TOWN 1/4 / TOW 2\ (mu:\^i'^ ( muj „ 



"mLo \ 1/4 



(113) 



as should be expected. 



3 Dissipative Systems 
3.1 Introduction 

In classical mechanics dissipation can often be adequately described by including 
a velocity dependent damping term into the equation of motion. Such a phe- 
nomenological approach is no longer possible in quantum mechanics where the 
Hamilton formalism implies energy conservation for time-independent Hamilto- 
nians. Then, a better understanding of the situation is necessary in order to 
arrive at an appropriate physical model. 

A damped pendulum may help us to understand the mechanism of dissipa- 
tion. The degree of freedom of interest, the elongation of the pendulum, under- 
goes a damped motion because it interacts with other degrees of freedom, the 
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molecules in the air surrounding the pendulum's mass. We may consider the 
pendulum and the air molecules as one large system which, if assumed to be 
isolated from further degrees of freedom, obeys energy conservation. The energy 
of the pendulum alone, however, will in general not be conserved. This single 
degree of freedom is therefore subject to dissipation arising from the coupling to 
other degrees of freedom. 

This insight will allow us in the following section to introduce a model for a 
system coupled to an environment and to demonstrate explicitly its dissipative 
nature. In particular, we will introduce the quantities needed for a description 
which focuses on the system degree of freedom. We are then in a position to 
return to the path integral formalism and to demonstrate how it may be em- 
ployed to study dissipative systems. Starting from the model of system and 
environment, the latter will be eliminated to obtain a reduced description for 
the system alone. This leaves us with an effective action which forms the basis 
of the path integral description of dissipation. 



3.2 Environment as Collection of Harmonic Oscillators 

A suitable model for dissipative quantum systems should both incorporate the 
idea of a coupling between system and environment and be amenable to an 
analytic treatment of the environmental coupling. These requirements are met 
by a model which nowadays is often referred to as Caldeira-Leggett model p^ , p3t 
even though it has been discussed in the literature under various names before for 
harmonic systems p^ , p5| , p6 27 and anharmonic systems psf . The Hamiltonian 



H = Hs + Hb + Hsb (114) 
consists of three contributions. The Hamiltonian of the system degree of freedom 

Hs^^ + V{q) (115) 

models a particle of mass m moving in a potential V. Here, we denote the 
coordinate by q to facilitate the distinction from the environmental coordinates 
Xn which we will introduce in a moment. Of course, the system degree of freedom 
does not have to be associated with a real particle but may be quite abstract. In 
fact, a substantial part of the calculations to be discussed in the following does 
not depend on the detailed form of the system Hamiltonian. 
The Hamiltonian of the environmental degrees of freedom 

N 



2m 

describes a collection of harmonic oscillators. While the properties of the envi- 
ronment may in some cases be chosen on the basis of a microscopic model, this 
does not have to be the case. Often, a phenomenological approach is sufficient 
as we will see below. As an example we mention an Ohmic resistor which as a 
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linear electric eleraent should be well described by a Hamiltonian of the form 



(116). On the other hand, the underlying mechanism leading to dissipation, e.g. 
in a resistor, may be much more complicated than that implied by the model of 
a collection of harmonic oscillators. 

The coupling defined by the Hamiltonian 



N N 



r.2 



n—1 n—1 

is bilinear in the position operators of system and environment. There are cases 
where the bilinear coupling is realistic, e.g. for an environment consisting of a 
linear electric circuit like the resistor just mentioned or for a dipolar coupling 
to electromagnetic field modes encountered in quantum optics. Within a more 
general scope, this Hamiltonian may be viewed as linearization of a nonlinear 
coupling in the limit of weak coupling to the environmental degrees of freedom. 
As was first pointed out by Caldeira and Leggett, an infinite number of degrees 
of freedom still allows for strong damping even if each environmental oscillator 

couples only weakly to the system [p2p^ . 

An environment consisting of harmonic oscillators as in ( |116| ) might be criti- 
cized. If the potential V{q) is harmonic, one may pass to normal coordinates and 
thus demonstrate that after some time a revival of the initial state will occur. 
For sufficiently many environmental oscillators, however, this so-called Poincare 
recurrence time tends to infinity p9[ |. Therefore, even with a linear environment 
irreversibility becomes possible at least for all practical purposes. 



The reader may have noticed that in the coupling Hamiltonian (117) a term 
is present which only contains an operator acting in the system Hilbcrt space 
but depends on the coupling constants c„. The physical reason for the inclusion 
of this term lies in a potential renormalization introduced by the first term in 



(117). This becomes clear if we consider the minimum of the Hamiltonian with 



respect to the system and environment coordinates. From the requirement 

- — = m„w„x„ - c„5 = (118) 

OXji 

we obtain ^ 

Xn = — • (119) 

Using this result to determine the minimum of the Hamiltonian with respect to 
the system coordinate we find 

dH dV ^ A cl dV 

oq oq mnOji, oq 

n—1 n—1 



The second term in (117) thus ensures that this minimum is determined by the 
bare potential V{q). 

After having specified the model, we now want to derive an effective de- 
scription of the system alone. It was first shown by Magalinskii |24| that the 
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elimination of the environmental degrees of freedom leads indeed to a damped 
equation of motion for the system coordinate. We perform the elimination within 
the Heiscnberg picture where the evolution of an operator A is determined by 



(121) 



From the Hamiltonian (114) we obtain the equations of motion for the environ- 
mental degrees of freedom 



Pn 



Pn 

m„ 



(122) 



and the system degree of freedom 



dV 
dq 



N 



N o 

c 



n=l 



P 

X = — . 

m 



(123) 



The trick for solving the environmental equations of motion (122) consists 
in treating the system coordinate q{t) as if it were a given function of time. The 
inhomogeneous differential equation then has the solution 



Xn{t) — a:„(0) cosiuJnt) + ^"^ ^ sin(a;„t) + 



ninUJn 



ds sin {t — s)) q{s) 



(124) 

Inserting this result into (|l23| ) one finds an effective equation of motion for the 
system coordinate 



mq - / ds V sin (a;„(t - s))q{s) + — + g V (125) 



N 



a;„(0) cos(aj„t) + ^"^^"^ sm{ujnt) 



By a partial integration of the second term on the left-hand side this equation 
of motion can be cast into its final form 



dV 

mq + m ds'-f{t — s)q{s) + —— — ^{t) 



(126) 



with the damping kernel 



N 



1 

= -5I^^cos(c^«0 

TO rUnUji 



(127) 
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and the operator-valued fluctuating force 



JV r / c \ fO) 

^(i) = c„ I x„(0) ^g(O) ) cos(tj„t) + — sin(cj„i) 

LV rrinUJ^ J m„w„ 



(128) 



The fluctuating force vanishes if averaged over a thermal density matrix of 
the environment including the coupling to the system 



m) 



B+SB 



TrB[^(t)exp(-/3(gB+gsB))] 
TrB[exp(-/3(i/B+i^SB))] 







(129) 



For weak coupling, one may want to split off the transient term m'j{t)q{0) which 
is of second order in the coupling and write the fluctuating force as [|o| 

C(t) - C(t) - rnj{t)q{0) . (130) 

The so defined force (^{t) vanishes if averaged over the environment alone 

TrB[C(i)exp(-/3i?B)] 



(CW)i 



TrB[exp(-/377B) 







(131) 



An important quantity to characterize the fluctuating force is the correla- 
tion function which again can be evaluated for ^ with respect to Hb + Hsb or 
equivalently for ( with respect to Hb alone. With (128) and (13C) we get the 
correlation function 



(C(i)C(0))B = V c„Q / f.T„(0) cos(a;„0 + sm{u;„t)] xi{0)) 

In thermal equilibrium the second moments are given by 



(a;„(0)x;(0))B = 6nl 



2mnUJn 



■ coth 



(P„(0)X((0))b = -^^&nl , 

so that the noise correlation function flnally becomes 



N 



(CWC(o))b = E 



2m„Cij„ 



coth 



V 2 



cos(a;„i) — isin(iu;„t) 



(132) 

(133) 
(134) 

(135) 



The imaginary part appearing here is a consequence of the fact that the o perators 
C(t) and C(0) in general do not commute. The correlation function ( 135 ) appears 
as integral kernel both in master equations as well as in the effective action 
derived below (cf. ( |l6^ ) and (p^). 
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It is remarkable that within a reduced description for the system alone all 
quantities characterizing the environment may be expressed in terms of the spec- 
tral density of bath oscillators 

N 2 

J(C0) = TT y ^ S{UJ - UJn) . (136) 

n— 1 

As an example, the damping kernel may be expressed in terms of this spectral 
density as 

7^ = - > — ^7 cos = — / ^ cos . 137 

For practical calculations, it is therefore unnecessary to specify all parameters 



mmi^n and c„ appearing in (116) and (117). It rather suffices to define the 
spectral density J{uj). 

The most frequently used spectral density 

J(tj) = mjuj (138) 

is associated with the so-called Ohmic damping. This term is sometimes em- 
ployed to indicate a proportionality to frequency merely at low frequencies in- 
stead of over the whole frequency range. In fa ct, in any realistic situation the 



spectral density will not increase like in ( 138 ) for arbitrar ily h igh frequencies. 
It is justified to use the term "Ohmic damping" even if (138) holds only be- 
low a certain frequency provided this frequency is much higher than the typical 
frequencies appearing in the system dynamics. 



From (137) one finds the damping kernel for Ohmic damping 

j{t) = 2j6{t) , (139) 



which renders ( |126| ) memory-free. We thus recover the velocity proportional 
damping term familiar from classical damped systems. It should be noted that 
the factor of two in ( |139D disappears upon integration in ( |126| ) since ( |137D implies 
that the delta function is symmetric around zero. 

At this point, we want to briefly elucidate the origin of the term "Ohmic 
damping" . Let us consider the electric circuit shown in Fig. ^ consisting of a 
resistance i?, a capacitance C and an inductance L. Summing up the voltages 
around the loop, one obtains as equation of motion for the charge Q on the 
capacitor 

+ (140) 

which shows that an Ohmic resistor leads indeed to memoryless damping. These 
considerations demonstrate that even without knowledge of the microscopic ori- 
gin of dissipation in a resistor, we may employ the Ohmic spectral density ( |138| ) 
to account for its dissipative nature. 
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R 



C 



L 



Fig. 9. LC oscillator with Ohmic damping due to a resistor R 



The spectral density (138) for Ohmic damping unfortunately diverges at high 
frequencies which, as already mentioned, cannot be the case in practice. Even 
in theoretical considerations this feature of strictly Ohmic damping may lead 
to divergencies and a cutoff is needed for regularization. One possibility is the 
Drude cutoff, where the spectral density 

J(^) = mjiu^^ (141) 

above frequencies of the order of wd is suppressed. The corresponding damping 
kernel reads 

7(i) =7^Dexp(-u;D|i|) . (142) 



This leads to memory effects in (126) for short times t < uj^^- For the long-time 



behaviour, however, only the Ohmic low frequency behaviour of the spectral 



density (141) is relevant. If a Drude cutoff is introduced for technical reasons, the 
cutoff frequency wd should be much larger than all other frequencies appearing 
in the problem in order to avoid spurious effects. 



The relation (136) between the spectral density and the "microscopic" pa- 
rameters implies that one may set Cn = mnUJn without loss of generality since 
the frequencies ujn and the oscillator strengths c^/2to„cl;„ can still be freely 
chosen. This special choice for the coupling constants has the advantage of a 
translationally invariant coupling |3l| ] 

n— 1 ^ ^ 

Furthermore, we now can determine the total mass of environmental oscillators 

= ^ Tdc.:^. (144) 

If the spectral density of bath oscillators at small frequencies takes the form 
J{lij) ~ w", the total mass of bath oscillators is infinite for a < 2. In particular, 
this includes the case of Ohmic damping where a free damped particle executes 
a diffusive motion. In contrast, for a > 2, the total mass is finite. In this case, 
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the particle will behave for long times like it were free albeit possessing a renor- 
malized mass due to the environmental coupling p^ . We emphasize that the 
divergence of the total mass for a < 2 is due to an infrared divergence and 
therefore independent of a high-frequency cutoff. 

It is also useful to express the potential renormalization introduced in (117) in 
terms of the spectral density of bath oscillators. From (136) it is straightforward 
to obtain 



N 



Si 

TT 



(145) 



This term is infinite for strictly Ohmic damping but becomes finite when a high- 
frequency cutoff is introduced. 



Finally, one finds for the noise correlation function (132) 



K{t) = (C(t)C(o)) 



B 



duj 



coth 



cos(ti)i) — isin(a;t) 



(146) 



In the classical limit, h 
expression 



0, this correlation function reduces to the real- valued 



K{t) = mkBT^it) , 



(147) 



where we have made use of (137). For Ohmic damping this implies delta corre- 
lated, i.e. white, noise. 

In the quantum case, the noise correlation function is complex and can be 
decomposed into its real and imaginary part 



K{t) = K'{t)+iK"{t) 



(148) 



Employing once more (137), one immediately finds that the imaginary part is 
related to the time derivative of the damping kernel by 



mSd7 
^ ' 2 dt 



(149) 



For Ohmic damping, the real part reads 



K'{t) 



(^^)%inh^ 



TTt 

hp 



(150) 



which implies that at zero temperature the noise is correlated even for long times. 
The noise correlation then only decays algebraically like 1/i^ much in contrast 
to the classical result (147). 



3.3 Effective Action 



In the previous section we had eliminated the environmental degrees of freedom 
to obtain the effective equation of motion (12(;) for the system degree of freedom 
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alone. This section will be devoted to a discussion of the corresponding procedure 
within the path integral formalism. 

We start to illustrate the basic idea by considering the time evolution of the 
full density matrix of system and environment 

W{q{,Xnf,qf,x'j^f,t) = J dqidq[dXnidx',^iK{qf,Xnf,t,qi,Xni,0) (151) 

X W{qi, Xni, q'i, x'„i, 0)K*{q'f, x'„f, t, q[, x^, 0) 

which is induced by the two propagators K. Here, the coordinates q and a;„ 
refer again to the system and bath degrees of freedom, respectively. The envi- 
ronment is assumed to be in thermal equilibrium described by the density matrix 
while the system may be in a nonequilibrium state p. If we neglect initial 
correlations between system and environment, i.e. if we switch on the coupling 
after preparation of the initial state, the initial density matrix may be written 
in factorized form 

W{qi,Xni,q'i,x'„i,0) ^ g()W|(x„i, <;) . (152) 

Since we are only interested in the dynamics of the system degree of freedom, 
we trace out the environment. Then the time evolution may be expressed as 



P('7f , q'l, ^) = y d(7idg( J((?f , q'l, t, gj, q[, 0)p{qi, ql) (153) 

with the propagating function 

J{q{,q'f,t,qi,ql,0) = J dXn{dXnidx'^iK{q{, Xnf,t, qi, Xni,0) (154) 

X W|(x„i, x'„i)K*{q'f, Xni, t, q{, x^, 0) . 

Here, the trace has been performed by setting Xni = x'^^ and integrating over 
these coordinates. The propagators may be expressed as real time path integrals 
while the equilibrium density matrix of the bath is given by a path integral in 
imaginary time. Performing the path integrals and the conventional integrals 
appearing in ( |154| ) one finds a functional depending on the system path. The 
important point is that this functional contains all information about the envi- 
ronment required to determine the system dynamics. 

For factorizing initial conditions, the propagating function J has been calcu- 



lated by Feynman and Vernon on the basis of the Hamiltonian (114). More 
general initial conditions taking into account correlations between system and 
environment may be considered as well p^ . 

Instead of deriving the propagating function we will demonstrate how to 
trace the environment out of the equilibrium density matrix of system plus envi- 
ronment. While this task is conceptually similar and leads to the same physical 
insight, it is considerably less tedious. 
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We start from the imaginary time path integral representation of the full 
equilibrium density matrix 



Wpiq.Xn.q'.x'^) = ^ / I I I T^^n ) cxp( -^5^[g,x„] ) (155) 

where the paths run from q(Q) = q' and a;„(0) — x'^ to q{hf3) = q and Xn{fiP) = 
Xn- The Euclidean action corresponding to the model Hamiltonian (114) reads 
in imaginary time 

5^ [g, x„] = Si [q] + Si [.T„] + [q, Xn] (156) 



with 



hp 



Si[q]^ dr(||2 + n9)) (157) 
Sl[xn]^ I ArY^'^iil+ulxl) (158) 



ji=i 

hl3 f N N 



5|b [q, Sn] = / dr -q- ^ c„5;„ + q^ ^ . (159) 

^0 V ^^^IrUnUJ,,! 



The reduced density matrix of the system is obtained by tracing over the 
environmental degrees of freedom 

Pl3{q,q') = ^rB{Wf3{q,Xn,q',x'j) 

I f _ f yJL /lp_ \ (160) 

= — / / _[J_dx„ * 11 IJcCn expl --5 [q,Xn]j 

where the circle on the second functional integral sign indicates that one has to 
integrate over closed paths a;„(0) = Xn{fiP) = Xn when performing the trace. 
The dependence on the environmental coupling may be made explicit by writing 



1 Jvqexp(^^^Si[q]^m (161) 



where the influence functional J-'[q] describes the influence of the environment 
on the system. Here, the partition function Z should not be confused with the 
partition function Zp of system plus environment. The relation between the two 
quantities will be discussed shortly. 

Since the bath oscillators are not coupled among each other, the influence 
functional may be decomposed into factors corresponding to the individual bath 
oscillators 



^ 1 

m = n ^^n[q] (162) 
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where 



1 



^" ~ 2sinh(ft/3c^„/2) ^^^^^ 

is the partition function of a single bath oscillator. The influence functional of a 
bath oscillator can be expressed as 



with the action 



hp 



dr- 



(164) 



(165) 



The partition function Z of the damped system is related to the full partition 
function Zp by the partition function of the environmental oscillators Z^ = 
Y\^^i Zn according to 2^ = Z^j Z-q. In the limit of v anish ing coupling, c„ = 0, 
the influence functional becomes = 1 so that (161) reduces to the path 
integral representation of the density matrix of an isolated system as it should. 
Apart from the potential renormalization term proportional to q^, the action 



(165) describes a driven harmonic oscillator. We may therefore make use of our 
results from Sect. 2_7. After analytic continuation t —ihf3 in ( p6| ) and setting 
Xi = X{ = Xn one finds for the classical Euclidean action 



:.E.ci 



m„UJn- 



1 



sinh(cj„r) + sinh(a;„ {hp - 
dr 



0) 



sinh{hl3LUn) 



Xnq{T) 



hp 
dr 

Tfln'-^n Jq Jo 

(-hp 



sinh(cj„(ft/3-r))sinh(tj„o-) 

do- . , — ^ q{vq{cr) 



2m 



n^n Jo 



(166) 



In view of the required integration over Xn one completes the square 



:.E.cl 



[q] 



1 



sinh(ft/3w„) 

„2 /.ft/5 



ft/3 



dr / daKn^T - a)q{T)q{a) 

10 



„2 i-np 

^ninUJi Jo 



(167) 



where Xn'^ does not need to be speci fied s ince it drops o ut a fter integration. 



The integral kernel appearing in ( 167 ) follows from ( 166 ) as 
2 cosh ra;„(^ - r) j 



sinh 



(168) 
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and therefore can be identified as the noise correlation function (135) in imagi- 
nary time 

Kn{T) = ■^(Cn(-iT)Cn(0))B - (169) 



The term in (|167| ) containing this kernel is quite unusual for an action. The 
double integral describes a nonlocal contribution where the system trajectory 
interacts with itself. This self-interaction is mediated by the environment as can 
be seen from the factor in (168). 

The integral kernel Kn{T) is only needed in an interval of length hp. Periodic 
continuation outside of this interval therefore allows us to expand the kernel into 
a Fourier series 

^n{T)^—T^ 2 r 2 exp(it^;r) 

1 — — 00 ^ 

2 oo 2^2 

eM^^ir) - TZZ-^ 1^ , ,2 I „2 exp(it^;r) 



c2 



J2 5{T-m-kn{T) , (170) 



where the Matsubara frequencies vi have been defined in (107). In (170), we split 
the kernel into two parts. The first term contains delta functions which lead to a 
local contribution to the action. Noting that due to the region of integration in 
(167) only half of the delta function contributes, this local term just cancels the 
potential renormalization ( |145| ). We are therefore left with the nonlocal kernel 

2 2 



fi(3mnUJ„ 



- — oo 



It can be shown that this kernel no longer contains a local contribution by writing 
dr / dakn{T - a)q{T)q{a) (172) 



hp 

dr / dafc„(r - a) {q{T) ~ q{<j)Y - [qirf + q[cTf) 



The first term is manifestly nonlocal because it contains the difference q{T)—q{a). 
Exploiting the symmetry of fc„(T), the second term can be expressed as 

1 fhp fT .hp .hp 

- dr / dafc„(T-a)(g(T)2+q(a)2) = / drqirf dafc„(a) . (173) 

^ Jo Jo Jo Jo 

Therefore, this term potentially could result in a local contribution. However, 
the time integral over the interval from to hj3 corresponds to the / = Fourier 
component which vanishes for A:„(r). As a consequence, the kernel knir) indeed 
gives rise to a purely nonlocal contribution to the action. 
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Now, we can carry out the Gaussian integral over x„ appearing in the influ- 
ence functional (164). With the action (167) we find 



I dr 



1 



2h 



hf) 



dcrfc„(T - a)q{T)q{(j) 



(174) 



The partition function Z„ arises from the fluctuation contribution and may be 
shown to be given by (163) for example by comparison with the uncoupled case 
c„ = 0. 

With (|162|) we finally obtain the influence functional 



T[q] = exp 



1 

'2h 



h(3 



dr 



dak{T — <7)q{T)q[(j) 



(175) 



with 



TV 



N 



V ' exp(iiy;T) 

^ ^ 4- 7/. 



2 



dw J(w) 



V ' 2 exp(ii^;r) 



(176) 



where we have made use of the spectral density of bath oscillators ( |13q ) to obtain 
the last line. 

The kernel hij) may be related to the damping kernel ^(t) by observing that 
the Laplace transform of the latter is given by 



7(^) 



r . , N / ^ 2 /-"^ dcj J{UJ) r . . ^ / ^ 

/ dicxp(— zf)7(t) = — / / di exp(— zf) cos(wt) 

Jo mjQ t: uj Jq 

duj J{w) z 



2 

m 



(177) 



In the first line we have employed the relation (137) between the damping kernel 
and the spectral density of bath oscillators. In view of ( |176D and (177) we can 
finally express the kernel as 



1 — — 00 



(178) 



For strictly Ohmic damping, this kernel is highly singular and we therefore in- 
troduce a Drude cutoff. The Laplace transform of the corresponding damping 
kernel is obtained from (142) as 



7(z) = 



7^D 
UJ-D + Z 



(179) 
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which reduces to 7(2;) = 7 for strictly Ohmic damping. Keeping the leading 
terms in the cutoff frequency cjd, the kernel now reads 



k{T)^m^u;u <5(r - nh^) - ^^^^ ^] + 0{uj^^) . (180) 



n=-oo sm 



hp 



For low temperatures, this gives rise to a long range interaction between different 
parts of the system trajectory. This reminds us of the algebraic decay of the real 



part (15C) of the noise correlation function K{t) and in fact it follows from our 



previous discussion that up to the periodic delta function appearing in (180) the 
kernel fc(T) equals —K{—\t). 

Summarizing this calculation we obtain the important result that the influ- 
ence of the environment on the system may be taken into account by adding a 
nonlocal contribution to the action. We then obtain the effective action 

-, php M 

S^M = S^[q] + -J^ dr da k{T - a)q{T)q{<j) (181) 



with fc(r) given by (178). The elimination of the environment within the real 
time path integral formalism, e.g. along the lines of the calculation by Feynman 
and Vernon p3] mentioned at the beginning of this section, would have led to an 



effective action of a structure similar to (181). An important difference consists 
in the fact that the propagation of a density matrix involves two paths instead 
of one. In addition, the integral kernel then of course appears in its real time 
version. 



4 Damped Harmonic Oscillator 



4.1 Partition Fiinction 

In this final section we will apply the results of the previous sections to the 
damped harmonic oscillator where exact results may be obtained analytically. 
The Hamiltonian describing system and environment is given by (114)-(117) 
with the potential 



Viq) 



2„2 



(182) 



As we have seen in Sect. B.2, there is no point in dealing with all the microscopic 
parameters present in the Hamiltonians (116) and ( |117|). Instead, it is sufficient 
to specify the spectral density of bath oscillators (136). In the following, we 
will mostly assume Ohmic damping, i.e. J{uj) = rwyu, and introduce a high- 
frequency cutoff of the Drude type (141) when necessary. 



In the general discussion of dissipative quantum systems we have concen- 
trated on imaginary time calculations and we will therefore try to take this 
approach as a starting point of the following considerations. Probably the most 
important quantity which can be obtained in imaginary time is the partition 
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function which in statistical mechanics can be viewed as a generating function 
for expectation values. Based on our previous discussion of the partition func- 



tion of the undamped harmonic oscillator in Sect. 2.9 and the effective action in 



imaginary time in Sect. 3.3, it is rather straightforward to obtain the partition 



function of the damped harmonic oscillator. 



According to (104) and (110) we can express the partition function of the 



undamped harmonic oscillator as 



n 



(183) 



We remind the reader that the denominator of the product stems from the 
fluctuation determinant associated with the Euclidean action 



hl3 



dr 



' 2 2 



(184) 



As we have seen in Sect. 3.3, the coupling to the environment leads to an addi- 



tional, nonlocal term to the action. Expanding the fluctuations in a Fourier series 
as we did on p. ^ and making use of the Fourier decomposition ( |178D of the 
integral kernel A;(t), we conclude that an additional term zy„7(i/„) appears in the 
fluctuation determinant. Modifying (183) accordingly, we find for the partition 
function of the damped harmonic oscillator 



h(3u}o 



n 



(185) 



For strictly Ohmic damping, we have 7(i^n) = 7. Since infinite products over 
terms of the form 1 -I- a/n for large n do not converge, we are forced to introduce 
a high-frequency cutoff in order to obtai n a finite result. One possibility is the 
Drude cutoff ( |l4^) with 7 given by ( |l79| ). 

In the following section we will try to extract some interesting information 
from the partition function and in the process will get an idea of where the 
difficulties for strictly Ohmic damping arise from. 



4.2 Ground State Energy and Density of States 

A thermodynamic quantity directly related to the partition function is the free 
energy which can be obtained from the former by means of 

F = -iln(Z). (186) 

In the limit of zero temperature, the free energy becomes the ground state energy 
of the undamped oscillator shifted due to the coupling to the environment. For 
the free energy, we find with ( |185|) 
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In the limit /3 — > cx) the spacing between the Matsubara frequencies Vn goes to 
zero and the sum turns into the ground state energy of the damped oscihator 
given by the integral 

It is particularly interesting to consider the case of weak coupling where 
connection can be made to results of perturbation theory. This will also help 
to understand the physical meaning of a ground state energy of a dissipative 
system derived from a free energy. An expansion of ( |l88| ) including terms of 
order 7 yields 

= A rd.lnfl + 4)+f rd.^^7M. (189) 

Evaluation of the first integral yields the expected result hu>o/2, i.e. the ground 
state energy of the undamped harmonic oscillator. The second integral repre- 
sents the shift due to the coupling to the environmental oscillators and can be 
expressed in terms of the spectral density of bath oscillators J(w). Recalling 



(177) one can perform the integral over v and the ground state energy (|189|) in 



the presence of damping becomes 

s,= ^J^ + J- r dcoJicu) , ^ , . (190) 
2 27rTO Jo ^ ^u;(wo + w) ^ ' 

In order to facilitate the physical interpretation, we rewrite this result as 

^0 = ^ - \ d^^(^) —. - - • (191) 



The first term of order 7 may be interpreted in analogy to the Lamb shift. 
There, an atomic level is shifted by creation and subsequent annihilation of a 
virtual photon as a consequence of the coupling to the electromagnetic vacuum. 
In our case, the atomic level is replaced by the ground state of the harmonic os- 
cillator and the environmental oscillators are completely equivalent to the modes 
of the electromagnetic field. The pictorial representation of this process is shown 
in Fig. ^ The coupling Hamiltonian (117) allows a transition from the ground 



state into the first excited state |1) by excitation of the j-th environmental os- 
cillator into its first excited state |lj). 

The energy shift associated with the diagram depicted in Fig. |l^ is given by 
second order perturbation theory as 

^o = fKilktoM)P (^92) 
^-^ —hujQ — huji 

where the denominator is determined by the energy of the intermediate state 
|1, Ij). With the matrix element 

(l,l,|ga;,|0,0) = - - (193) 
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Fig. 10. The ground state energy of the harmonic oscillator is shifted by a transition 
to the first excited state accompanied by a virtual excitation of the j-th environmental 
mode 



and the relation (|136|) for the spectral density of bath oscillators we get 



An=—, 



2'KmuJo Jq 



dujJ{iLj] 



(194) 



which is just the first term of or der 7 in (191). As we know, the bilinear coupling 
Hamiltonian appearing in (192) gives rise to a renormalization of the potential 
which has been taken care of by the second term in the Hamiltonian (117). 
The result (194) contains this potential renormalization since only the bilinear 
coupling term has been considered. In (191), which results from the full Hamil- 
tonian, this effect is subtracted off by the second term under the integral in (191) 
as can be verified by comparison with (145). 

It is obvious that for strictly Ohmic damping with J{uj) = 777.70; the correc- 
tion (194) and with it the ground state energy (191) will display an ultraviolet 
divergence which is due to the unphysical behaviour of the spec tral density J(uj) 
at large frequencies. Assuming a Drude cutoff we find with (141) to leading order 
in the cutoff frequency lud the finite result 



Aujo 27r V ^0 / 



(195) 



The negative first term corresponds to the potential renormalization which is 
no longer present in the ground state energy Sq- The second term, on the other 
hand, is positive and thus leads to an increase of the ground state energy. 

Not only the ground state energy can be derived from the partition function 
but one may also formally introduce a density of states p{E) of the damped 
system according to [p5| 



/•OO 

Z(/3) = / dEp{E)exp{^l3E) 
Jo 



(196) 



Inversion of the Laplace transformation allows to determine p{E) from the par- 
tition function according to 



p{E) 



1 
27ri 



d/3Z(/3)exp(/3£:) 



(197) 
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where the constant c has to be chosen such that the hnc of integration is to the 
right of all poles of Z{f3). 

Once a high-frequency cutoff for the spectral density of bath oscillators is 



specified, the inverse Laplace transform in (197) may be evaluated either numer- 
ically or by contour integration. The second approach leads to a series which 
again has to be evaluated numerically |35| . However, it is not necessary to intro- 
duce a cutoff provided we shift the energy by the ground state energy Eq which 
in fact is the only divergent quantity in this problem. Such a shift may be per- 
formed by considering Zexp(/3eo) instead of Z itself. To demonstrate that this 
procedure renders the cutoff irrelevant, we will restrict ourselves to the limit of 
weak damping and large cutoff considered before even though a more general 
treatment is feasible. 

In a first step we decompose the infinite product appearing in the partition 



function (185) with Drude cutoff (179) into a factor where the limit uju oo 



can safely be taken and a factor still containing the cutoff frequency 

n— 1 ^ n— 1 ^ ' / , \ 

It is the last product which has to be analyzed with care because it vanishes in 
the limit cj/j — > oo. To leading order in 7 one finds 

/ 00 \ oc 

7^D \ \ - 71^0 



\n=l 



Un{v,, + Wd) ) VniVr, + Wd) ^^^^^^ 

where we have introduced the digamma function ||3^ 

00 

^■(1 + ^) = -C -f y . (200) 

' ^ n(n + z) ^ ' 

Here, C ~ 0.577 ... is the Euler constant. With the leading asymptotic behaviour 
+ z) ^ ln(z) for large arguments z, the partition function for large cutoff 
frequency becomes 

Z - ^ ft (^] "''"'^ (201) 



On the other hand, we know from our result (195) that the ground state energy 
diverges to leading order with (?i7/27r) ln(a;D)- Therefore, multiplication of the 
partition function with exp(/3eo) leads indeed to an expression with a finite value 
in the limit of infinite cutoff frequency. 

Before taking a look at numerical results, we remark that the partition func- 
tion contains a pole at /3 = with residue 1 / hujQ. This pole represents the Laplace 
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6- 




5 10 15 20 

{E - eo)/hujo 



Fig. 11. The density of states ( |l97| ) of the damped harmonic oscillator is shown for 
7/2ajo ~ 0.05. A delta function contribution at E — eo has been omitted. The dashed 
line marks the average density of states 1/hujo 



transform of the constant 1 / Hluq and therefore is related to the average density 
of states which takes the value of the undamped case where the energy spacing 
between adjacent levels is hajQ. In Fig. ^ we present the density of states for 
weak damping, 'y/2u!o = 0.05. A delta function contribution at E — Eq has been 
omitted. Due to the weak damping we find well defined peaks which are close 
to the energies expected for an undamped oscillator. With increasing energy 
the levels become broader. For stronger damping, only the lowest levels can be 
resolved and a level shift induced by the damping becomes visible. 

The behaviour of the level widths shown in Fig. ^ is consistent with the 
result of a perturbative treatment. According to Fermi's golden rule, the width 
of the n-th level is given by 

27r °° 

-f^n = T2 X! !("- + 0)| 6{~uJo-u;j) 

^ 0=1 (202) 

I 1 2 

+ \ {n — l,lj\cjqXj\n,Q)\ 5{ijJo-ujj) 

where we have already taken into account that the matrix element of the dipole- 
type coupling only connects the state n to its nearest neighbours. Because of 
energy conservation the first part of the sum never contributes. With the matrix 
elements 

(n-l,l,|qx,|n,0) = — (203) 

2,[mmjU}t)UJj) ' 

we thus find for the width of the n-th level in terms of the spectral density of 
bath oscillators 

71 

r„ = J{uq) . (204) 

For Ohmic damping, J{lo) = rwyuj, we finally get 



(205) 
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As observed before, the levels broaden with increasing damping strength 7 and 
level number n. We remark that one can demonstrate by a semiclassical analysis 
of other one-dimensional potentials that it is indeed the level number and not 
the energy which is decisive for the level width . 



4.3 Position Autocorrelation Function 



In the introduction we have mentioned fluctuations as one of the effects arising 
from the coupling to an environment. Even if a system is i n the rmal equilibrium 
with its environment, fluctuations due to the noise term (128) will be present. 
The appropriate quantity to describe this phenomenon are equilibrium correla- 
tion functions like the position autocorrelation function 



Cit) = {qit)q{0)) = Triq{t)q{0)pp) 



(206) 



From this quantity one can derive all other equilibrium correlation functions of 
the damped harmonic oscillator as is discussed in Q. 

We now want to determine this correlation function by first calculating its 
imaginary time version and start with the Euclidean action 



hp 



dr 



2 2 
-u^q 



h / ATF{T)q{T) . 

Jq 



hp i-hp 
dr / Aak{T — (j)q{T)q{a) 

(207) 



The second term accounts for the coupling to the environment as we have shown 
in Sect. 3.3, In addition, we have included the third term corresponding to an 
external force in imaginary time. This constitutes a useful trick which will allow 
us to determine the correlation function by variation with respect to this force 



{q{T)q{a)) = fi^Tv 



SF{t) SF{a 



■P 



(208) 



As we know already from Sect. 2.7, the force does not appear in the fluctua- 
tion part. It is therefore sufficient, to restrict our attention to the classical path. 
The classical equation of motion following from variation of the action (207) 



mq{T) 



hp 



dcrk{T — o-)q{a) — mujQq{T) = F{t) 



(209) 



is most conveniently solved by Fourier transformation on the interval from to 
hf3. Introducing the Fourier transforms 



n— — 00 



qn exp(ii^„T) 



and 



F{r) 



m 



^ /nexp(ii^„r) 



(210) 



(211) 
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and making use of the Fourier representation (178) of A;(r) for Ohmic damping, 
wc find for the Fourier coefficients of the classical solution 



fn 



(212) 



Inserting this result into the Fourier representation of the action (207) 



S^[q] = ^ [('^^ + ^'^"1 + ^ohnl-n + fnq-n + f-nQn] 



(213) 



yields the classical Euclidean action 



m 



E 



fn f— r 



(214) 



or equivalently 



2mhl3 



E 



1 



/ dr / daF{T)F{a) exp (iiy„(T - a)) . 
Jo Jo 



(215) 



Since the external force appears only through this action in the exponent of the 
equilibrium density matrix, we may easily evaluate the functional derivatives 
according to (208). For the position autocorrelation function in imaginary time 
we thus find 

C{r) = ^ E ^"Tf 2 (216) 

Unfortunately, the real time correlation function cannot be obtained simply by 
replacing the imaginary time r by it where i is a real time. For negative times t, 
the sum in (216) would not converge. We therefore have to perform an analytic 
continuation to real times in a more elaborate way. 

The idea is to express the sum ( ^16[ ) as a contour integral in the complex 
frequency plane. To this end, we need a function with poles at frequencies co = 
ii/„, n = ~oo, . . . ,oo with a residuum of one. This requirement is satisfied by 
hf3/[l — exp(— /t/^Lj)]. By integration along the contour shown in Fig. |l^(a) we 
find 



dw- 



exp(— wr) 



c+ 



1 — exp{—hPio) —uj'^ 



-1-I-2-E 



exp(ii/„r) 



(217) 




Fig. 12. The analytic continuation of tlie imaginary time correlation function 
(g(r)g(O)) to real times mak es use of the integration contours depicted in (a) and 
(b) to obtain (217) and (218), respectively 



Similarly, an integration along the contour shown in Fig. |T2|(b) leads to 



, hp exp(— wr) 

doj- 



exp(— fi,/3a;) —uj^ — ijuj + ujq 



. TT , „ . exp(iiy„T) 
1^ + zvri > 

^0 ^n+7i^n + 



(218) 



Subtracting (218) from (217), the imaginary time correlation function can be 
expressed as 

1 >^ exp(iz^„T) 



E 



' n— — oo ' ^ /oi n\ 

h f°° "fuj exp(— wr) 

' duj- 



''^'^ J-oo ('^^ ~ ^o)^ + 7^ijJ^ 1 — cxp(— fi,/3a;) 

Now we may pass to real times by the replacement r — > it to obtain the real 
time correlation function 



^(*) = dc. , . S ../"'^7't^^ • (220) 

f^'-'"' J-oo ("^ ~ '^o) + J 1 — exp(— apw) 



Physical insight into this result can be gained by considering the Fourier 
transform of this correlation function 

/oo 
dtexp(iujt)C{t) (221) 
-oo 

which may be related to the dynamic susceptibility x(<^) of the damped harmonic 
oscillator. According to the Ehrenfest theorem the equation of motion for the 
expectation value of the position agrees with the corresponding classical equation 
of motion. In the presence of an external force F{t), the latter reads 

mq + ui'-fq + mujlq = F{t) . (222) 
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The solution of this equation may be expressed in terms of the response function 
Xit) as ^ 

qit) = f dsx{t~s)Fis) (223) 

<J —OO 

which by means of a Fourier transformation becomes 

Q» = x{^)F{lu) . (224) 



With the equation of motion (222) the dynamic susceptibiUty of the damped 



harmonic osciUator is then found to read 



X(^) = - 2 _u 2 • (225) 



Together with (|220| ) and ( p2l| ) one finally obtains the relation 

1 — exp(— ft/3w) 

which represents an example of the so-called fluctuation-dissipation theorem 
. Here, the position autocorrelation function describes the fluctuations while 
the imaginary part of the dynamic susceptibility x" can be shown to determine 
the energy d issip ated under the influence of the external driving F{t). While 



the relation (226) is exact for linear systems like the damped harmonic oscilla- 
tor considered here, it still holds for nonlinear systems within linear response 
theory. There, the response to the external force is considered in leading order 

perturbation theory . 

It is instructive to consider the real time correlation function ( ^20| ) in more 
detail. We first decompose the correlation function C{t) into its real and imagi- 
nary part, or equivalently, into its symmetric and antisymmetric part 



with 



and 



C{t) ^ S{t) + iA{t) , (227) 



S{t)^\{{q{t)qm + {q{Q)q{t))) (228) 



A{t)^~'-{{q{t)qm~{q{Q)q{t))) . (229) 



From mm we find 



S{t) r dc.'-^ ^ cothf ^ ) cos(c^t) (230) 

271771 J _oo (OJ"^ — Wq) + 7 



and 



h f°° 

A{t) = —- / d.; 
27rm 



—z TjTT, siTi{ujt) . (231) 
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Apart from Planck's constant appearing in the prefactor, the antisymmetric part 
is purely classical. In fact, one demonstrates within linear response theory the 
general relation 

Xit) = ^{[q{t),q{0)mt) = -lA{t)@{t) . (232) 

Therefore, our statement follows as a consequence of the Ehrenfest theorem 
which ensures that the response function x{t) of the damped harmonic oscillator 
is classical. 

More interesting is the symmetric part ( |230| ) of the correlation function C{t). 
There exist two different types of time scales determined by the poles of the 
integrand in ( ^3Cl| ). One set of four poles at frequencies 

uj = ± { oj zt i 



((i±i|) (233) 



corresponds to the characteristic frequencies of a damped harmonic oscillator 
with the oscillation frequency 



^=[^0-^) (234) 

shifted by the damping. In addition, there exists an infinite number of poles at 
imaginary frequencies ivm n — — cx3, . . . , oo depending on the temperature of the 



heat bath via the Matsub ara f requencies defined in ( [1071 ). With these poles, one 



can evaluate the integral ( 23C ) by means of a contour integration to obtain 
/z_ ui /o^ [sinh(^/3tZ>) cos(tZ>0 + sm{hp-//2) sin(c^|^|)] 

_ ^ t^nexp(-^-n|t|) (2^^^ 

'"^ ^1 + ^o)^ - l^'^n ' 

The sum in the second line becomes important at low temperatures k^T <^ 
h'y /At:. In order to discuss this quantum effect, we focus in the following discus- 
sion on the case of zero temperature. Then, the poles of the hyperbolic cotangent 
in ( ^30|) become dense and form a cut on the imaginary axis which has conse- 
quences for the long-time behaviour of th e correlation function. In the limit 
/3 — > oo the sum in the second line of ( |235| ) turns into an integral. Noting that 
the long-time behaviour is determined by the behaviour of the integrand at small 
arguments we find 

Sit) r dx^^^^t^ - — ^4 . (236) 



Instead of the usual exponential decay we thus find an algebraic decay. 

In the limit of vanishing damping, the imaginary part of the dynamic suscep- 
tibility appearing in the integrand of (23C ) turns into a sum of two delta functions 
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located at loq and —loq. For weak but finite damping, the delta functions broaden 
into Lorentzians 



'JLU 7/1 1 



Lj2)2_|_^2^2 4<I> - + 72/4 (cj + + 72/4 



(237) 



corresponding to the four poles ( 233 ) and one can assume that only the integrand 
in the neighbourhood of these poles is relevant. Within the so-called Markov ap- 
proximation, one then replaces the contributions to the integrand which depend 
only weakly on frequency by their values at Q. As we will see, in contrast to 



(236) the correlation function S{t) at zero temperature then no longer decays 



algebraically. It is interesting to discuss the reason for this discrepancy. 



To this end we go back to the integral representation ( 230 ) of the correlation 
function S{t). In a first step, we apply the so-called rotating wave approximation 
(RWA) which consists in neglecting the Lorentzian located at — w, i.e. the second 



term in (237) ||4^. Physically speaking, we neglect processes where the system 
is excited into an energetically higher state by loosing energy to the driving 
force or where the system decays to a lower state by absorbing energy. For finite 
temperatures, we now have 

SRWAit) = ^ _ [ duj- — ^- — 277 C0th( ^ ) cos(wt) . (238) 

STrmw (01-^)^-1-7 /4 V 2 / 

Within the Markov approximation, we replace the hyperbolic cotangent by its 
value a.t to ~ uj. In the limit of zero temperature this leads to 

S'RWA,Markov(i) = T T / duj- —5— 577 COs(wt) 

Snmuj [uj - uY + Y/^ (-239) 

^ t- \ / 7 



■ cos{ijjt) exp^- 



-t 



Amui V 2 

We thus find an oscillation with a frequency shifted due to the environmental 
coupling and an exponential decay in contrast to the algebraic decay ( |236| ). 

This difference can be understood by taking a closer look at the Markov ap- 
proximation. In Fig. [13 the Lorentzian and the hyperbolic cotangent appearing 



in (238) are schematically shown as full and dashed line, respectively. In order 
to obtain ( p39D , we have approximated the hyperbolic cotangent by a constant. 
However, in doing this, we have replaced an antisymmetric function by a sym- 
metric one which can only yield a non- vanishing result together with the rotating 
wave approximation made above. As a consequence, the area shaded in grey has 
been taken with the wrong sign. The idea was that this difference should be 
small and arising from frequencies far from ui. However, due to the large extent 
of a Lorentzian, it is of order 7 and, in addition, it replaces an exponential decay 
by an algebraic decay. 



At zero temperature the difference between (238) and (23£) becomes 

^ = 'S'RWA(i) — '5'rWA, Markov (i) 

h f° 7 (240) 

= -2- / doj- —5 :r—COs(ujt) 

8nmJ_^ (a;-cj)2+72/4 
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Fig. 13. The full and dashed lines represent the L orent zian and hyperbolic cotangent, 
respectively, which contribute to the integrand in (238) 



which may be expressed in terms of integral sine and integral cosine functions. 
For our purpose it is sufficient to note that for long times, the difference indeed 
decays algebraically as 

We remark that the factor of two relative to the result ( |236| ) arises because no 
rotating wave approximation has been made in deriving the latter. 

In the previous discussion, we have been concerned with an effect of order 
7 which in the spirit of the Markov approximation should be considered as 
small. However, the Markov approximation changes the long-time behaviour of 
the equilibrium correlation function S{t) qualitatively and the absence of an 
algebraic decay, even with a small prefactor, may be relevant. For truly weak 
coupling, the damping constant 7 should be the smallest frequency scale. Apart 
from the usual weak coupling condition 7 <C t'.'o we also have to require 7 ^ 
k^T/h which at low temperatures may become difficult to satisfy. 
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